Small scale rotational disorder observed in epitaxial graphene on SiC(OOOl) 
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Interest in the use of graphene in electronic devices has motivated an explosion in the study of this remarkable 
material. The simple, linear Dirac cone band structure offers a unique possibility to investigate its finer details by 
angle-resolved photoelectron spectroscopy (ARPES). Indeed, ARPES has been performed on graphene grown 
on metal substrates but electronic applications require an insulating substrate. Epitaxial graphene grown by 
the thermal decomposition of silicon carbide (SiC) is an ideal candidate for this due to the large scale, uniform 
graphene layers produced. The experimental spectral function of epitaxial graphene on SiC has been extensively 
studied. However, until now the cause of an anisotropy in the spectral width of the Fermi surface has not been 
determined. In the current work we show, by comparison of the spectral function to a semi-empirical model, 
that the anisotropy is due to small scale rotational disorder (~ ± 0. 15°) of graphene domains in graphene grown 
on SiC(OOOl) samples. In addition to the direct benefit in the understanding of graphene's electronic structure 
this work suggests a mechanism to explain similar variations in related ARPES data. 



Interest in the single layer of hexagonally coordinated car- 
bon atoms, known as graphene, has been intense ever since the 
discovery of its unusual electronic properties An under- 
standing of the electronic properties is essential if graphene 
applications are to be realised. Angle Resolved Photoemis- 
sion Spectroscopy (ARPES) measurements provide the most 
direct method to investigate the electronic band structure, 
however this technique requires large well defined samples. 
ARPES have been performed on graphene on metal sub- 
strates, however electronic applications will require insulat- 
ing or semiconducting substrates. A significant number of 
ARPES studies have therefore been performed on graphene 
grown epitaxially on SiC by thermal decomposition yl [3[], a 
preparation method which is a very promising candidate for 
applications because the production process can be scaled up, 
uses common semiconductor processing steps and provides 
large scale uniform layers. 

This material can be considered an important testbed for ex- 
citing solid state physics, and interest has turned to fine details 
of its band structure, as observed by ARPES. In this context 
we compare ARPES measurements to a semi-empirical pho- 
toemission model. Comparison of the model to experimental 
data accounted for almost all features of the electronic struc- 
ture, with the exception of an anisotropy in the spectral width. 
We show that this anisotropy is explained by a ± 0.15° ro- 
tational disorder of the graphene domains. Fig. [U(a) shows 
the well known experimental spectral function obtained from 
epitaxial graphene on SiC(OOOl), ~ -0.5 eV, which was 
demonstrated to sit on top of the so-called buffer layer with 
(6V3 x 6%/3)i?30° periodicity (hence forth 6^) Hi- The 
anistropy in the Fermi surface intensity Jsj], the increase in the 
intensity between the Fermi surface and ~ 200 meV and 
the offset of the bands above and below the Dirac crossing, 
Ed , J3-E1] have all been successfully described. Until now 



the anisotropy in the spectral width of the Fermi surface seen 
in Fig. [TJ(a.i) as an increased width in the vertical, k y (KK), 
direction when compared to that in the horizontal, k x (TK), 
direction has not been explained. 

This variation is not confined to the Fermi surface, but ap- 
pears across all energies. This is important for a detailed anal- 
ysis of the spectral function in many contexts such as gap 
formation, the examination of many-body interactions and 
the investigation of the the Berry's phase in graphene. The 
current work provides a simple explanation of this feature, 
in terms of small scale rotational disorder, which is essen- 
tial to the discussion of the electronic structure of graphene. 
In fact many papers have used intensity and line-width vari- 
ations in ARPES data to discuss the existence of various 
quasipa rticles||9- 151 and in discussions of the Berry's phase 
111 l[ll6lll7ll in graphene. The current work therefore provides 
an invaluable insight into ARPES spectral variations and in 
the understanding of these more exotic features. 



ARPES were obtained at the Maestro end station (SES- 
R4000 analyzer) at beamline 7 of the Advanced Light Source, 
Lawrence Berkeley National Laboratory. Spectra where taken 
at ~ 20 K, at < 2xl0 -10 Torr using 95 eV photons giving 
an overall resolution of ~ 25 meV and 0.01 A -1 . Epitaxial 
graphene on SiC(0001) was prepared either by annealing in 
ultra-high vacuum (UHV)|2|] (vacuum grown) or by anneal- 
ing in Ar (argon grown) using a custom build reactor| 18 1. 
The decomposition of SiC in Ar leads to an improved crys- 
talline quality yfl. 



The bare band used in the semi-empirical ARPES model is 
a first nearest neighbour tight binding (FNN TB) model fit to 
experimental graphene on 6\/3 data. This models the n bands 
via the relations: 
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with a lattice constant, a = 2.46 A, and the fitted parameters, 
7o=-3.24 eV and So=0.0425 eV, are those found by Bostwick 
et al. 11711 . The final fitted parameter, E2 P , is the offset of the 
Dirac energy, Ed , from the Fermi level due to doping of the 
graphene by the substrate and is determined by comparison to 
the ARPES measurements. Broadening of the bare band is in- 
troduced to the model through the self energy via the spectral 
function relation fljj 
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The self energy is determined using the semi-empirical 
method of Bostwick et al. J20I1 . where the linewidth of the 
ARPES momentum density curves (MDCs) are used to de- 
termine the imaginary component of the self energy, which 
is Hilbert transformed to get the real component. This ex- 
perimental self energy is then used to recreate the experimen- 
tal data and a self-consistent fitting used to further refine the 
self energy. An additional Gaussian broadening term is added 
to the current model to account for experimental broadening, 
dE = 25meV. 

Rotational disorder of graphene on metal substrates such 
as Cu(lll)[21] and Pt(lll)0| has been shown to lead to 
significant anisotropic spectral broadening, therefore we at- 
tempted to model the much smaller anisotropic broadening 
seen in Graphene on SiC with a similar small rotational dis- 
order. To model the effect of rotational disorder several 3D 
spectral functions, rotated around the V point by small angles, 
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FIG. 1. Experimental (a) and semi- empirical model (b) Fermi sur- 
faces (i) and spectral functions in the TK (ii) and KK (iii) directions 
of epitaxial graphene on SiC(OOOl). 



are summed together. This is illustrated, for a 2D Fermi sur- 
face, in Fig. |2](b). The non-rotated Fermi surface is the solid 
arc around the K point, two further Fermi surfaces (dashed 
arcs) are shown rotated by a small amount around T. The ef- 
fect on the Fermi surface is that it appears broader in the KK 
direction than in the TK direction. Examples of the Fermi sur- 
face intensity profiles in the KK and TK directions for three 
different regions on a sample are shown in Fig. |2](c). Fitted 
Fermi surface data presented in Fig. [3]and Fig. |4]are obtained 
by taking similar intensity line profiles of the Fermi surface in 
1 degree steps around the K point with the TK direction being 
set as the zero angle. The Lorentzian width is set to the width 
of the Fermi surface in the TK direction (0. 15 A -1 ), ensuring 
the Gaussian line width describes the variation in width of the 
Fermi surface. 

In order to investigate the anisotropy in the width of the 
Fermi surface, experimental spectral functions at a number 
of K points in the graphene Brillouin zone ( see Fig. [2] (a)) 
were obtained. Fitted Fermi surface data were then extracted 
from these spectral functions and are presented in Fig. [3] and 
Fig. |4] Spectral functions obtained using the semi-empirical 
model described in the methods section were also analysed in 
a similar manner and are overlayed on the data in Fig. [3] and 
Fig. ED 

The experimental data obtained from the first Brillouin 
zone (position A, B and C from Fig. |2]l is presented in Fig. 
[3] This data is compared to a semi- empirical model which in- 
volves 15 rotational domains equally spaced between ±0.15°. 
The size (peak position) and intensity (peak amplitude) of the 
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FIG. 2. Diagram showing the (a) Brillouin zone boundaries for 
graphene with the 1st, 2nd, 3rd and 4th Brillouin zone K point radial 
lines (dashed, blue, lines) of graphene, (b) schematic of a selected 
K point showing the effect of small scale rotational disorder and (c) 
selected experimental and semi- empirical model intensity profiles in 
the KK and TK directions. In (a) the Fermi surfaces indicate the K 
points at which measurements where undertaken. 
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FIG. 3. Fitted Fermi surface data for vacuum grown (solid light, orange, lines), argon grown (solid dark, green, lines) and semi- emperical 
model (dashed, blue, lines) for the three 1st Brillouin K points indicated in Fig. [2] Position(i), amplitude(ii) and Gaussian Linewidths (iii) 
obtained from Lorentzian-Gaussian fits to the Fermi surface are presented. In all cases the Lorentzian width is set to 0.15 A -1 ensuring that 
the Gaussian width describes the variation in the width of the Fermi surface. 



Fermi surface is well described by the model. The anisotropy 
in the width of the Fermi surface (variation in the peak Gaus- 
sian linewidth) is clearly evident by the "peanut" shape in Fig. 
[5] (a.iii), (b.iii) and (c.iii) for both vacuum grown and argon 
grown samples. The "peanut" shape is only produced in the 
semi-empirical model by considering rotational disorder. The 
size of the "peanut" lobes is directly proportional to the ro- 
tational variation of the domains, with ±0.15° providing the 
best agreement. 

Due to the six-fold symmetry the first Brillouin zone K 
points are all equivalent with respect to the proposed rota- 
tional disorder. This is not the case for the higher order Bril- 
louin zones (BZ2, BZ3 and BZ4 in Fig. |2]i where the axis 
of rotation no longer lies in the TK direction. Investigation 
of the higher Brillouin zone experimental spectral functions 
(presented in Fig. |4| then provides an important test of the ro- 
tational disorder model. The size (peak position) and intensity 



(peak amplitude) of the Fermi surfaces at each of the higher 
Brillouin zone locations are again well described by the rota- 
tional disorder model (15 domains, ±0.15° rotational spread). 
In particular Fig. |4](b.ii) indicates that the experimental Fermi 
surface intensity has a maximum at ~ 190° (dash dot dot, light 
green, line) rather than the expected 180°, which is predicted 
by the rotational disorder model (dash dot, red, line). 

The "peanut" shape of the Fermi surface anisotropy (Fig. [4] 
(a.iii) ) is observed in the second Brillouin zone, However this 
becomes more oval shaped in the third (b.iii) and fourth (c.iii) 
Brillouin zones as the axis of rotational disorder no longer 
aligns with the intensity anisotropy axis J^. This variation is 
also well described by the rotational disorder model. 

Further evidence for the rotational disorder is found in ex- 
perimental Fermi surface intensity profiles (dashed lines) in 
the KK and the TK directions (Fig. [2] (c) taken from three 
distinct regions on an argon grown sample. The intensity pro- 
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FIG. 4. Fitted Fermi surface data for vacuum grown (solid light, orange, lines), argon grown (solid dark, green, lines) and semi emperical 
model (dashed, blue, lines) for the higher order Brillouin zone K points indicated in Fig. [2] Position (i), amplitude (ii) and Gaussian Linewidths 
(iii) obtained from Lorentzian-Gaussian fits to the Fermi surface are presented. In all cases the Lorentzian width is set to 0.15 A" 1 ensuring 
that the Gaussian width describes the variation in the width of the Fermi surface. The red and green lines in b ii indicate the position of the 
maximum intensity. 



file from the semi- empirical model, with the correct rotational 
disorder, is overlayed in red (solid line). The spectra were ob- 
tained from the first Brillouin zone K point and show a similar 
width and shape in the TK direction (Fig. |2] (cii), (civ) and 
(cvi)). In contrast three distinct line shapes are observed in 
the KK direction (Fig. |2](c.i), (c.iii) and (c.v)). 

The most commonly observed lineshape across several 
samples is shown in Fig. 12 (ci) and is well described by the 
rotational disorder model (15 domains within the ±0.15° ro- 
tational spread). The other two spectra are rare compared to 
the rotationally disordered spectra but provide a significant in- 
sight. In Fig. |2] (c.iii) an asymmetric line shape is observed, 
which is described by a model including only 2 rotational do- 
mains (±0.15°) with an intensity ratio of 2:1 between them. 
The final region (Fig. [2] (c.v)) has a similar lineshape in the 
KK and TK direction, corresponding to a region containing 



only a single rotation. 

Comparison of experimental Angle Resolved Photoemis- 
sion Spectra from "vacuum grown" and "argon grown" epi- 
taxial graphene on SiC(0001) to a semi-empirical model con- 
firms the existence of a < ±0.15° rotational disorder of the 
Fermi surface. This disorder is attributed to a number of ro- 
tated graphene domains within the 50/im photon beam size. 
Experimental data from regions of a vacuum grown sam- 
ple which show only 1 and 2 rotational domains are also 
presented, however these are rare compared to the rotation- 
ally disordered spectra. Importantly the scale of the rotation 
(±0.15°) is much smaller than is possible to determine from 
modern imaging techniques, like LEEM/PEEM and STM, and 
the common diffraction techniques, LEED and RHEED. It is 
therefore shown that detailed analysis of ARPES features can 
give information on small scale structure variations not possi- 
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ble until now. 
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